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ABSTRACT 


The deformation of an unstable baroclinic wave as a result of mean 
horizontal wind shear is studied with the use of a five-layer primitive 
_ equation model. Two basic wind profiles were tested: 1) uniform zonal 
flow, 2) jet type zonal flow. The jet profile produced a cold front 
and a warm front, with the warm front being the more intense. The 
uniform zonal flow produced a cold front. The cold fronts were pro- 
duced by horizontal shearing deformation fields and the warm front was 
produced by horizontal stretching deformation fields. The formation of 
the warm front was also enhanced by strong horizontal convergence. The 


occlusion process was not observed in the experiments. 
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I. INTRODUCTION 


The concept of atmospheric fronts has been accepted for many years, 
yet at this time there is no mathematical theory which can explain the 
formation of a front. Williams [1967] has indicated that the basic 
problem is finding a process which can produce a small-scale motion 
field from one which initially contains only large-scale motions. 

A number of studies have been conducted evaluating mechanisms 
which are contributing factors in the formation of a front from an 
fefitial large-scale baroclinic wave. Stretching deformation wind 
fields have —_ studied by Stone [1966] and Williams and Plotkin 
[1968], while shearing deformation wind fields have been studied by 
Arakawa [1962] and Williams [1967]. Stone [1969] indicates that these 
effects should be of about the same importance. Arakawa [1962], 
Edelmann [1963], Williams [1965, 1967] and $kland [1969] have studied 
atmospheric frontogenesis with vertical deformation fields with the 
general result that a realistic front is produced within a reasonable 
time period. Bushby and Timpson [1967] and Bushby [1969] have studied 
the Apianies of fronts and frontal rainfall with a primitive equation 
model utilizing real data with good results in a twenty-four hour 
forecast of rainfall. 

The purpose of this study is to investigate the deformation of a 
unstable baroclinic wave as a result of horizontal shear in the mean 
wind field. The primitive equation model is basically the model used 
by Harrison [1969] with a slightly modified lower boundary condition 


and different geographical and theoretical areas of application. 
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Il. THE PRIMITIVE EQUATIONS 


Consider the primitive equations in the following form, with 
heating, friction and curvature terms neglected and with the coriolis 


parameter a linear function of y: 


- =- Lt) - os fv, (1) 
& = - L(v) - 2 sofia. (2) 
= - 10), (3) 
w_ ou ov 
Ca X(2y 
Op ae Op GO : a 
where 
L(s) =< CULE RE tee ae) + = (ws ), 
and S=u,v, 9. 
At p = Pp. 
: = -S ws) -L wh +9 S+S -o @. (6) 


Equations (1) through (3) and equation (6) are prognostic 
equations. Equations (1) and (2) are the momentum equations, 
equation (3) is the first law of thermodynamics and equation (6) 
is the lower boundary condition. The diagnostic equations (4) and 
(5) are the mass continuity and hydrostatic equations, respectively. 

Equation (6), the lower boundary condition was obtained by 


setting df/dt = gw = 0, which implies that w = 0 at the lowest level. 
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Phillips [1963] has shown that the error resulting from the neglect of 


vertical motion at the lowest level is very small. 
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Itt. FINITE DIFFERENCES 


Equations (1) through (6) were solved numerically by introducing 
finite differences in x, y, p and t. The space finite differences 
employed were a variation of a scheme constructed by Arakawa [1966] 
[also see Lilly, 1965]. The basic scheme is designed to conserve 
total energy and also to conserve the squares of u, v and 9 with 
respect to the advection terms. 

The vertical domain is divided into five layers with the arrange- 


ment of variables as shown in Fig. l. 


mb Variable Level (k) 


0 ——— we0 —— 10 

—— 100 ——— u,v,f,e ——— 9 —— 
200 3———— ~~ 8 

nO — —_- —_ u,v,9,9 — — — 7 ——— 
400 —_——_— ———_—__--—— 6 

—— 500 —-——-— u,v,f,e@ ——-— 5 —— 
600 9———___—— a e 

—-—— 700 = —~—-_ a, ,9 — — — 32 oe 
agg SS SS ee 

—— 900 —— — —_s vf, 0 — — ad 
1000 + ———— u,v,f,¢,0 ——————__ 0 


Fig. 1. Distribution of variables in the vertical. 


The finite difference form of the L(s) operator can be written 
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The linear space derivatives are of the general form 


“F- 
ah 
The hydrostatic equation (6) can be written in finite difference 


form following Lorenz [1960] as 
g = zi Pag ig an WF EDe 
ihe kane LAK 2 


PY’ Plcb oN 
“p — me | 
The form of vertical finite differences as specified in the L(s) 
Bereat or and the Bees eabic equation were modified at level zero in 
the following manner. 
1. @, u and v were linearly extrapolated from levels one and 
three to level zero. 
2. 0of/dp at level zero was approximated by a one-sided 
| difference. 
Centered time differences were used throughout this study with the 
exception of the first time step which is forward. To prevent sepa- 


ration of solutions, the program was restarted every twelve hours with 


a forward time step. 
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The following basic steps were used in integrating the given 
system of equations. 

1. The potential temperature and wind fields at the odd- 
numbered levels and the geopotential of the 1000 mb surface are 
prescribed initially. 

2. The mass continuity equation is integrated in space from 
the highest level to the surface to obtain w. 

3. The hydrostatic equation is integrated in space from the 
surface to the highest level to obtain @. 

4. The momentum equations are integrated in time to obtain 
u and™v ate time t + At. 

5. The first law of thermodynamics is integrated in time to 
jeeatn @ at time t + At. 

6. The equation for the lower boundary condition is integrated 
in time to obtain the 1000 mb geopotential at time t + At. 

7. Steps 2 through 6 are repeated for the number of time steps 


desired. 
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IV. BOUNDARY CONDITIONS 


The boundary conditions are specified in a manner that conserves 
mass and energy within the grid domain. At the northern and southern 
boundaries the normal component of velocity is zero and periodicity is 
enforced at the eastern and western boundaries. Fig. 2 illustrates the 


horizontal grid domain. 


i=l i=M 
ies spud 2. mune N. boundary —— ee 52 
eie@ é > ¢ 6 e © e a 
lee © @ e ® oe 
' ’ 
W. boundary E. boundary 
oa ) ¢ é a 
: : 
Fn id a 2 os °@ oo» @ ° — = 
a = 2 s ne e ae a s we 
i=1 S. boundary i= 
jel 


j=l 
Fig. 2. Horizontal grid | 

The east-west boundary condition forces cyclic continuity on all 
parameters. This is accomplished by requiring the values in the two 
columns on the right to be the same as the values in the two columns 
on the left. This condition is specified by 

“1,j  M-1, 3 "2,4. -M, 3? 
where 
s = u,v,$,9. 

The north-south boundary condition allows flow along the boundary 
but does not allow flux across the boundary. This is accomplished by 
imposing the following condition on the v term of the advection operator, 
ve + Vv. = 0 7 2 = 0 


claret 2 i,N i,N-1 
or 


lags 


The values of geopotential and potential temperature in the outer row 
are set equal to the value in the next interior row. This is speci- 
fied by 

| S. ie se Sy =s, (7) 


where 


As a result of v being set equal to zero between the two outer grid 
points along the northern and southern boundaries, values of u and w 
are not required outside the boundaries. 

The upper boundary condition of w = 0 eliminates flux through the 
top of the domain. The lower boundary condition has been specified in 
equation (6). 


Various boundary conditions were tested (see Appendix A) before 


acceptance of the above method. 
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V. INITIAL CONDITIONS 


The initial fields of u, v and 9 were specified in order to 
generate a baroclinic zonal current. A small amplitude disturbance 
in the geopotential field; which is independent of pressure, was 
superimposed on the basic zonal flow. The mean zonal current (u) has 
the following forms; 

1) jet profile 
of 5) int (®) 
Op O W 


2) uniform zonal profile 


where Ou/dp is constant. The initial data fields were generated in 
the following steps. 


1. The geopotential at p = p, was calculated from 


” L. 2T1X ey (®) 
p , Ex U I COS ee sin \y (8) 
2. The potential temperature fields were specified from 


p n 
P a Se 5 
/ du Se (, . 
Q R Sp Jiasts ae ur Sp (p Pads 


Ns 
where 0@/dp is a constant. 


For uniform zonal flow 


7 2 
Ji= fy + BLY ioe 


ie, 


For jet flow 


I= fy + By*/2 
a (2m 
oF EF: GT aa W 


g y 2 
_lofy. . Ww ary) _ a gj ary | 
2 E pea cos (40 mT sin (0 


S. The hydrostatic equation was then solved for the geopotential 
at all.required levels. 

4, The initial wind fields were then calculated from the geo- 
potential fields using the geostrophic approximation. Thus, the 
initial fields were in geostrophic balance at time zero. 

The basic current is zero at the lowest level and increases at a 
rate of 5 mps/100 mb to a maximum of 45 mps at 100 mb. Therefore 
du/ dz decreases with increasing altitude which reduces the baroclinic 
instability in the upper layers. This approximates the effects aie 
- tropopause on the dynamics. 

The initial disturbance amplitude was 10 meters or .125 mb. The 
small initial amplitude was chosen to prevent the excitation of large 
amplitude gravity waves that occur since the initial wind field was 
geostrophic. 

The value of 80/dp used was -.376°C/mb. This value is approxi- 
mately one-half the value determined by Gates [1960] for winter over 
the United States, but it has the advantage of forcing the fronts to 
develop more rapidly and thus to reduce the computer time required. 
The value used approximated the standard atmosphere quite well in the 
lower levels, see Table I, but is less stable than the standard atmos- 


phere in the upper levels. The stratosphere is clearly not represented 
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by this temperature field, however the dynamical presence of the 


stratosphere is represented by u. 


Height (meters) 


Pressure (mb) Standard 
Jet Atmosphere 





Table I. Average initial height of constant pressure surface. 


The beta plane approximation was made, where the coriolis parameter 
is a linear function of y. The quantities f and 8 were computed at 
O 
45° north latitude. 


Various combinations of grid size, time step and wind profile were 


tested. Table II is a summation of these various experiments. 
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VI. RESULTS 


All experiments were conducted using the values of ANWR Gp 58/dp 
and initial disturbance amplitude as specified by equation (8) in the 
initial condition section of this paper. Table II gives a listing of 
the specific variations in each of the following experiments. 

Experiment No. 1. This is the basic jet flow experiment with a 
Prid sine of 200 km in x and y. The initial 1000 mb height analysis 
is given in Fig. 3 and the initial height and potential temperature 
analysis for the 900 mb and 700 mb surfaces are given in Figs. 4 and 
5. The growth of the systems from day one through day four at the 900 
mb surface is given in Figs. 6 through 9. By the end of three days a 
cold front and a warm front have formed and the temperature gradients 
are large by the fourth day. Note that the warm front contains larger 
gradients than the cold front. Large convergence (not shown) is associ- 
ated with the warm front. The 700 mb analysis of height and potential 
temperature for days three and four is given in Figs. 10 and 1l. The 
cold and warm fronts have the proper tilt in the vertical, with a 
slope between 1/125 and 1/250, which is computed at day three and for 
both fronts. The growth of the surface pressure centers, divergence 
and vorticity at 900 mb is given in Table III. It should be noted in 
Table III the asymmetry between the maximum and minimum values of each 
quantity. 

Experiment No. 2. This is the basic uniform zonal flow experiment 


with a grid size of 200 km in x and y. The initial 1000 mb height 


analysis is given by Fig. 3 and the initial height and potential 


a8) 





Table III. Experiment No. l. 


temperature fields for the 900 mb and /00 mb surfaces are given in 

Figs. 12 and 13. The growth of the systems from day one through day 
four at the 900 mb surface is given in Figs. 14 through 17. By day 
three a cold front has formed to the west of an apparent warm front. 


However there is a single vorticity maximum in the region and there 


= 


is no evidence of a warm front at higher levels. The front appears 
to be an occlusion but it did not evolve from two separate fronts with 


associated vorticity maxima. It was also noted that in the northern 


region of the low there was a packing of potential temperature lines 
which somewhat resembles the warm front in the jet flow experiment. 
The 700 mb pressure analysis of height and potential temperature for 
days three, and four are given in Figs. 18 and 19. The front tilts to 
the west in the vertical with a slope between 1/60 and 1/120 at day 
three. The growth of the surface pressure centers, divergence and 
vorticity at 900 mb is given in Table IV. 

Experiment No. 3. This experiment is the same as experiment l 
except that the grid length was reduced to 100 km in x and y. The 


results are basically the same as in experiment 1, with the exception 
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1000 mb Divergence Vopr icadtby 


-6 = 
DAY centers (m x LO x aO 


Low | High Max Min Max Min 


fe fe [ea [oe [ie [oe 


| 4 -550 340 7 VR, -43.7 ow) “4.5 


Table IV. Experiment No. 2. 

that the growth of the systems is now slightly faster with a slightly 
larger vorticity. The 900 mb height and potential temperature analysis 
for days three and four are given in Figs. 20 and 21. Note that the 
temperature gradients are larger in this experiment at four days than 
“in experiment 1. This indicates that the larger truncation error with 
the larger grid size inhibits the energy cascade to small scales. 

Experiment No. 4. This experiment is the same as experiment 2 
except that the grid size was reduced to 100 km in x and y. A front 
similar to the one in experiment 2 developed, but in addition, a 
gravity wave with a wavelength of approximately 188 km was super- 
imposed on the height analysis. The gravity wave was initially 
observed at day two close to the north and south boundaries, and by 
day three had completely covered the domain. It was apparently excited 
by the large shear at the boundary as generated by equation (7). 

Experiment No. 5. This experiment is the same as experiment l 
except that Ay was reduced to 100 km and W was reduced to 2000 km. 


Cold and warm fronts developed within three days that were similar to 
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the ones generated in experiment |. Superimposed on the long wave 
pattern was a very small amplitude gravity wave with a wavelength of 
approximately 188 km. 

Experiment No. 6. This experiment is the same as experiment 2 
except that Ay was reduced to 100 km and W was reduced to 2000 km. 

The results are the same as in experiment 2, but in addition, a 
gravity wave of small amplitude and approximate wavelength of 188 km 
was superimposed on the long wave pattern. 

Experiment No. 7. This experiment is the same as experiment l 
except that Ay was increased to 400 km and W was increased to 8000 km. 
See Figs. 22 and 23 for day three 900 mb height and potential tempera- 
ture analysis. By day three a cold and warm front, similar to the 
results of experiment 1, had formed with the disturbance being limited 
to the central region of the domain. 

Experiment No. 8. This experiment is the same as experiment 2 
except that Ay was increased to 400 km and W was increased to 8000 km. 
" Gravity wave propagation, with an approximate wavelength of 400 kn, 
“from the north and south boundaries was too strong to permit an 
evaluation of the long wave system. 

In none of the experiments was there an indication of a front at 
500 mb or above. All energy calculations (see Appendix B) indicated 
that the potential energy was three orders of magnitude greater than 
the kinetic energy and that there was a small energy sink in the 


system up to four days. 
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Minimum contour 100 meters, contour interval 
10 meters. 


Fig. 3. Initial 1000 mb height analysis. 
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Minimum contour 800 meters, interval 30 meters. 


Fig. 22. Experiment No. 7, 72 Hour 900 mb height analysis. 
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Minimum Isotherm 270°A, interval 5°C. 


Fig. 23. Experiment No. 7, 72 Hour 900 mb potential temperature 
analysis. 
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VII. CONCLUSIONS 


The time evolution of a baroclinicly unstable disturbance was 
studied for two basic zonal currents. With horizontal wind shear, 
cold and warm fronts developed and are similar to the results of 
Edelmann [1963] and Qkland [1969]. With no horizontal wind shear, 

a cold front formed which was similar to the results of Williams 
[1969]. This cold front resembled an occluded front at the surface 
but did not have the life history of an occlusion. In general, the 
warm front formed from the stretching deformation discussed by Stone 
[1966] while the cold fronts formed from the shearing deformation 
discussed by Williams [1967]. In the jet case, the at front was 
further enhanced by strong convergence through the action of the 
vertical deformation field. The lack of frontogenesis at upper levels 
is explained by the reduction of the vertical wind shear at those 
levels. The classical occlusion process was not observed in any of 
the experiments. 

Further studies with this model should be conducted with the 
addition of vertical and horizontal diffusion and with the addition 
of latent heat of condensation. Other basic currents should be 


tested, including those with an explicit tropopause. 
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APPENDIX A 


EXPERIMENTS WITH BOUNDARY CONDITIONS 

Before acceptance of the north-south boundary conditions as 
specified in Chapter V, various other schemes were tested with 
limited success. A sample of these schemes with the results follows. 

1. The boundary was placed at the outer grid points with the 
values of all variables determined from a linear extrapolation from 
the next two interior grid points. The integration became unstable 
after approximately 120 six-minute time steps. 

2. The boundary was again placed at the outer grid points and a 
special one-sided difference modification was made to the basic diag- 
nostic and prognostic equations. Strong unrealistic wind shear 
developed along the boundary and the numerical integrations eventually 
became unstable. 

3. Cyclic boundary conditions were applied in the north-south 
as well as in the east-west directions. The discontinuity in the 
geopotential field forced the use of a complete wave in both the 
north-south and east-west directions. This doubled the number of grid 
points in the domain and made the size of the computer program too 


large for practical use. 
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APPENDIX B 


ENERGY COMPUTATIONS 
In pressure coordinates, the total potential energy may be 


written as 


Po 


I CT dp da, Beall 


Ao 


[= 

tg 

TT 
co |e 


where "dA" represents an area increment. Now, substituting the 
definition of potential temperature in (B.1), the expression for 


the total potential energy becomes, 


mB! oe) dp dA. B.2 


In pressure coordinates, the total kinetic energy is written 


as 


Es 
phon aie striata on dp dA. B.3 
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ee 
> —, 
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